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I.  INTRODUCTION 


Bimetallic  disks  with  a small  initial  curvature  have  been  used  as  the 
temperature- sensing  and  control  element  of  thermostatic  switches  for  several 
decades.  The  first  detailed  analysis  of  the  thermoelastic  stability  of  such  a 
disk  was  presented  in  1953  in  a two-part  paper*'  in  Part  1,  Wittrick* 
derived  a general  theory  of  the  thermoelastic  deformation  of  such  a disk  and 
investigated  the  condition  for  its  snap-through  instability;  in  Part  II.  Witt  rick, 

Myers  and  Blunden^  obtained  numerical  solutions  to  the  nonlinear  boundary- 

value  problem  derived  in  Part  I with  the  aid  of  a differential  analyzer  (also 

known  as  an  analog  computer).  Wittrick's  analysis  provided  the  groundwork 

for  understanding  the  operation  of  the  bimetallic-disk  (BMD)  thermostat  by 

showing  that,  if  the  initial  central  height  of  the  disk  is  greater  than  a critical  V 

value— approximately  twice  the  disk  thickness— > the  disk  will  respond  to 

temperature  changes  with  very  positive  snap-through  actions  at  two  fixed 

temperatures,  one  on  heating  and  the  other  on  cooling.  This  positive  snap- 

through  action  is  a desirable  feature  in  tripping  electric  contacts  since  it 

minimizes  arcing  and  thus  prolongs  the  service  life  of  the  contacts. 

In  recent  years,  the  BMD  thermostat  has  attracted  particularly  favorable 
interest  in  the  aerospace  industry.  The  benefit  of  its  application  in  a space- 
craft or  satellite  system  is  obvious:  the  present-day  BMD  thermostat  can  be 
remarkably  lightweight  and  small  in  size*,  sealed  for  protection  from  con- 
tamination. and  yet  sturdy  in  construction  for  proper  functioning  in  a shock- 
and- vibration  environment.  Furthermore,  it  operates  solely  on  the  thermo- 
elastic principles  and  requires  no  external  power  or  a gravitational  field  to 

*A  commonly  used  model  weighs  only  8.5  g (0.3  os),  measures  1.  5 cm  (0.6  in.) 
in  diameter  and  length. 

1 Witt  rick,  W.H.,  "Stability  of  a Bimetallic  Disk,  Pert  L " Quart.  J.  M.ch. 

Appl.  Math..  VI  <19531.  p.  IS. 

2Wlttrlck.  W.H.,  DM.  Mycrt,  and  W.  R.  Pluad.n,  "Stability  of  a Bimetallic 
Dlak,  Part  11,  ■ Quart.  J.  Mach.  Appl.  Math..  VI  U953),  p.  26. 


trip  the  switch  (as  required,  respectively,  by  a solid-state  temperature 
sensor  or  a mercury -bottle -type  thermostat).  All  these  features  can  qualify 
the  miniaturized  BMD  thermostat  as  uniquely  suitable  for  space  applications 
if  its  design,  manufacture  and  testing  can  be  controlled  to  ensure  the  high 
level  of  reliability  required  in  space  programs.  Unfortunately,  the  current 
state-of-the-art  for  both  theoretical  understanding  and  practical  fabrication 
of  the  BMD  thermostat  has  not  yet  reached  such  an  advanced  stage,  partly 
due  to  the  fact  that  since  the  publication  of  Witt  rick's  work*  nearly  a quarter 
of  a century  ago,  very  little  analytical  effort  has  been  made  to  further  the 
understanding  of  the  BMD  thermostat  operations,  ^s  evidenced  by  the  scarcity 
of  applicable  references  in  the  technical  literature. 

It  has  been  well-known  among  both  manufacturers  and  users  of  the  BMD 
thermostat  that  many  such  units  can  exhibit  anomalous  switching  behaviors. 

One  such  undesirable  switching  anomaly  is  a drift  in  both  upper  and  lower 
switching  temperatures.  In  1972,  Henry  and  Coffin*  presented  a combined 
analytical  and  experimental  study  addressing  this  drift  problem.  Based  on 
Wittrick's  analysis,  extensive  numerical  calculations  were  performed  to  show 
that  high  cyclic  and  mean  stresses  may  develop  in  a typical  bimetallic  disk 
during  its  snapping  cycles.  They  also  compared  their  analytical  results  with 
strain  gage  measurements,  obtaining  reasonably  good  agreement.  These 
results  stron*iy  support  their  proposition  that  high  cyclic  and  mean  stresses 
in  the  disk  are  the  major  contributing  factor  of  the  drift  in  itt»  snapping  tem- 
peratures during  prolonged  cyclic  operation.  As  a result  of  their  parametric 
study,  a design  procedure  was  proposed  in  Ref.  ? to  limit  the  maximum 
stresses  in  the  disk,  thus  alleviating  the  drift  problem. 

Recent  experience  of  extensive  applications  of  the  miniaturised  BMD 
thermostat  in  space  systems  has  brought  to  light  some  new  problems  associated 
with  its  switching  behavior.  The  most  frequently  encountered  switching 
anomaly  is  the  so-called  temperature  "creepage"  or  "dithering"  phenomenon 
(a  thermostat  exhibiting  such  an  anomaly  is  sometimes  called  a "creeper" 

^Henry,  M.  F.,  end  L,  F.  Coffin,  Jr.,  "An  Investigation  of  Switching  Stresses 
in  Bimetal  Disks,  " Int.  J.  Mech.  Sci.,  M.H972),  p.  345. 
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or  a "ditherer").  These  terms  have  been  used  inte rchangeably*  to  describe 
either  of  the  following  two  conditions:  (1)  the  failure  of  an  assembled  unit  to 
respond  to  temperature  changes  with  immediate  positive  snap  action  of  the 
disk  at  one  of  the  set-temperatures  during  manufacturers'  or  users'  accep- 
tance test,  or  (Z)  a sudden  deviation  of  the  switching  temperatures  of  a unit  in 
service  from  its  original  set-temperatures,  resulting  in  a very  narrow 
switching  band.  The  former  condition  usually  results  in  high  rejection  rates 
during  acceptance  tests  and  thus  in  increased  manufacturing  cost.  The  latter 
case  can  be  far  more  Consequential,  since  such  events  are  usually  charac- 
terised by  either  a hesitant  contact  or  a series  of  frequent  openings  and 
closings  of  the  contacts  and  thus  may  induce  excessive  arcing  or  stress  cycling: 
these  effects  often  result  in  switch  malfunction  or  shortened  service  life. 

In  the  following  sections.  Wittrick's  analysis*  will  be  first  extended  to 
include  the  effects  of  a central  force,  representing  the  interaction  between 
the  disk  and  an  elastic  switching  armature;  then  numerical  results  will  be 
presented  to  show  that  the  type  of  switching  anomalies  discussed  above  may  be 
attributed  to  the  improper  regulation  of  such  an  interaction  force. 


In  this  report,  the  two  terms  "crcepage"  and  "dithering"  will  be  used  with 
slightly  different  meanings  instead  of  interchangeably  (see  Section  V - B for 
detailed  discussions). 


II.  BIMETA  LLIC  - DISK  THERMOSTAT  DESIGN  AND  OPERATION 


There  ere  many  possible  designs  of  e BMD  thermostat,  and  it  can  be 
used  to  operate  either  a heat  source  or  a heat  sink  to  provide  necessary 
thermal  protection.  For  simplicity,  it  will  be  assumed  that  all  the  thermo- 
stats discussed  herein  are  used  to  control  an  electric  heater  so  that  it  is  on 
when  the  contacts  are  closed  and  off  when  separated. 

The  basic  designs  of  the  BMD  thermostat  may  be  categorised  into  two 
general  types,  (a)  and  (b),  as  illustrated  in  Fig.  1.  In  Fig.  1(a),  the  contacts 
are  closed  with  a built-in  preload  between  them  when  the  disk  is  disengaged 
from  the  " strike r-pln,  " while  in  Fig.  1(b)  the  contact  force  is  supplied  by  the 
disk-pin  interaction.  For  applications  in  a space  system,  the  type  (a)  design 
is  somewhat  preferable  to  type  (b)  for  the  following  reasons: 

1.  Since  the  contact  force  in  type  (a)  is  preselected  during  design 
and  may  be  calibrated  during  assembly,  while  that  in  type  (b) 
is  solely  determined  by  the  disk  ac't<  1,  the  latter  inevitably 
has  a higher  random  variance  fro n *ne  switch  to  another  due 
to  manufacturing  limits  of  dimensional  precision  and  material 
uniformity.  This  random  variance  in  the  contact  force  of 
type  (b)  switches  may  result  in  undesirable  variance  in  their 
electrical  characteristics  and  service  life. 

2,  The  elastic  cantilevered  armature  in  the  type  (a)  design  has  an 
additional  firm  support  at  all  times  during  Us  operation,  i.e., 
either  the  contacts  or  the  striker-pin  (or  both)  are  engaged, 
while  that  in  type  (b)  has  no  additional  firm  support  when  the 
switch  is  in  the  off  position*  Therefore,  when  vibrations  or 
dynamic  transients  are  present,  the  armature  In  type  (b)  wiki 
be  more  susceptible  to  vibratory  excitation  that  may  cause 
accidental  contact  closing  or  arcing. 

For  these  reasons,  the  following  discussion  and  analysis  will  primarily 
concern  the  type  (a)  design;  however,  their  extension  to  type  (b)  is  straight- 
forward. 
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Fig.  1.  Schematic  Drawing  of  Bimetallic -Dish  Thermostat 


At  this  point,  consider  the  operation  of  a typical  BMD  thermostat  as 
shown  in  Fig.  1(a).  For  the  configuration  shown,  the  switch  is  on  and  the  disk 
is  in  a "prebuckling " equilibrium  state  (i.e.,  on  the  same  side  as  its  stress- 
free  state,  such  as  points  A or  G in  Fig.  2),  with  the  high-expansivity  metal 
layer  on  the  concave  side.  As  the  disk  temperature  rises  due  to  either  heater 
output  or  ambient  temperature  rise,  the  disk  will  deform  until  it  reaches  an 
unstable  position1  (point  B in  Fig.  2)  at  the  upper  snap  temperature  0^.  At 
this  temperature,  the  disk  snaps  through  toward  a "postbuckling"  equilibrium 
state  (point  C in  Fig.  2).  However,  during  the  actual  thermostat  operation, 
the  disk  cannot  reach  the  equilibrium  state  Casa  free  disk;  it  will  first  hit 
the  striker-pin  at  point  C^,  deform  the  armature  until  the  contacts  separate 
at  point  and  then  settle  to  an  equilibrium  state  represented  by  point  C'. 
Note  that  in  order  to  determine  this  new  equilibrium  state  C\  the  disk-pin 
interaction  force  must  be  included  in  the  analysis.  However,  since  this  point 
represents  the  position  of  the  system  at  rest,  it  suffices  to  treat  this  inter- 
action as  a static  force  (i.e.,  the  impact  and  dynamic  transient  components  of 
the  force  may  be  ignored). 

Further  examination  of  the  nature  of  the  armature  deformation  reveals 
that,  if  the  disk  and  armature  materials  remain  in  the  linear  elastic  range 
during  switch  operations,  the  disk-pin  interaction  force  will  be  a bilinear 
function  of  the  disk  center  height  h,  as  shown  in  the  upper  part  of  Fig.  2.  The 
stiffness  of  the  armature  at  the  striker-pin  will  assume  a higher  value  kj 
between  points  Cj  and  than  the  value  k^  for  points  beyond  :ue  to  the 
additional  support  at  the  contact  point.  If  the  armature  is  very  rigid  with  a 
soft  rotational  spring  at  the  root,  the  bilinear  function  approaches  a step 
function  in  the  limit  (kj/k2-*oo),  and  points  Cj  and  C2  coincide.  If  the  arma- 
ture is  a uniform  cantilever  beam  with  the  striker-pin  at  the  middle  point, 
k | / * 32/7,  For  the  actual  sample  switches  analyzed,  the  armature  is  a 
thin  plate  with  a specially  designed  shape  (sometimes  plated  for  better  con- 
ductivity), and  the  ratio  k | /^2  ran8ea  between  five  and  seven. 
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m.  ANALYSIS 


louaider  the  bimetallic  dl»k  shown  in  Dig , L Kxeept  tor  the  intmduc 
tion  of  a central  toYce  V and  the  necea »a rv  edge  support,  the  prot>lem  ia 
identical  to  that  treated  in  Kef.  I;  the  additional  edge  aupport  determine*  an 
i nt  eg  rat  i on  constant  in  the  deflection  \v  which  differ*  from  that  in  Ket.  I only 
hv  a rigid  body  uiaplacement , It  . The  notation  lined  herein  it*  caHentially 
the  same  a*  that  m Ket,  l with  minor  modification*  to  conform  to  current 
prevailing  notation  lined  in  ahell  theory  luce  Symbol*  aectionl, 

A.  liOVKKNINr.  KUUA TIONS  KOK  A DISK  WITH  A 
i KNTKAI.  KOH ( K 

Since  the  ha»ic  kinematic  relation*  and  conatitut  ive  equation*  for  the 
problem  under  conaide  ration  are  the  name  a*  Kq*.  (11  through  (1*11  of  Ket.  I 
and  are  alan  available  in  Ket.  they  will  not  he  repeated  here.  The  oipia 
tion*  of  equilibrium  including  the  central  force  effect  are^ 
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where  w i*  the  total  di»tance  from  a point  on  the  reference  *urtace  ot  the  di*k 
to  the  bane  plane  p.e,,  vs  sv  ^ t w^,  where  w ^ *)»ecifie*  the  initial  ahape  aiul 
\v  ^ denote*  the  the  rmoe  la  *t  tc  detlectionl,  and  a po*itive  V'  i*  in  the  ponitive 
g direction.  Note  that  Kq.  (.*»  i*  the  reault  of  integrating  the  differential 
equation  ot  equi librium  in  the  t.  direction,  with  the  central  force  F appearing 
m the  integration  con*tant. 
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KeiHMiier,  I*!,,  “Symmetric  bending  ot  Shallow  Sliell*  ot  Kevolution,  " 

,1.  Math.  U Merit.,  7 (1^8).  p.  1 .*  I . 
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Upon  elimination  of  the  shear  stress  resultant  Qr>  Eqa.  (2)  and  (3) 
combine  to  give 
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Equation  (I)  is  the  same  as  Kq.  (14)  of  Kef.  1,  while  Eq.  (4)  differs  from 
Eq.  (15)  thereof  only  by  the  additional  term  E/27T. 

Denote  the  elastic  rotation  of  a normal  to  the  disk  as  a new  variable/?; 
then,  within  the  scope  of  a thin  shell  theory 
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where  8 l dw  /dr.  The  equivalent  flexural  rigidity  l>*  and  the  equivalent 
ro  o 

Poisson's  ratios  V.  and  are  now  defined  as  follows; 
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where  the  quantities  B,  r,  t),  D\  l’*,  ancl  D'  are  the  same  as  in  Kef.  1.  The 
same  elimination  procedure  used  in  Kef.  1 reduces  the  basic  set  of  equations 
to  the  following  two  coupled  differential  equations  |cf.,  Eqs.  (In)  and  (171  of 
Kef.  I); 
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The  two  boundary  conditions  at  the  disk  edge  are  the  same  as  that  in  Ref.  1; 
i.e., 


and 


N =0  at  r = a 
*• 


+ V, 


D r 


at  r = a 


(8) 


(9) 


where  & is  the  disk  temperature  measured  from  a reference  temperature  at 
its  stress-free  state,  and 


Q*  * Q 


PC 

B 


in  which  Q and  P are  also  defined  in  Ref.  1.  The  other  two  boundary  condi- 
tions at  the  disk  center  must  be  modified  to  account  for  the  singularity  intro- 
duced by  the  central  force  F.  If  F is  assumed  to  be  small,  so  that  it  causes 
only  a linear  perturbation  of  the  thermoelastic  deformation,  the  linear  singular 
solution  for  such  a shallow  shell,  similar  to  those  given  in  Refs.  S through  8, 
applies: 


liTo  Nr  ■ SR*  “fc  - V C r"o  ,n(s)  + «"«*•  <‘0) 


!‘Zt  ■sfc,',ro  ["”(;)]  <“> 

Reiaaner,  E,,  "Streaaea  and  Small  Diaplacementa  of  Shallow  Spherical  Sheila, 
Part  I.  " .1.  Math.  I.  Phya..  25  (19471.  p.  80. 

^Relaaner,  E.,  "St  re  a tea  and  Small  Diaplacementa  of  Shallow  Spherical  Sheila, 
Part  11.  w J.  Math.  «.  Phya..  25  (1947).  p.  279. 

^Fliigge,  W.,  and  R.  Etting,  "Singular  Solutiona  for  Shallow  Sheila,  * Int,  J. 
Solida  L Struc..  8 (1972).  p.  227. 
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Simmonds.  J.G.,  end  M.  R.  Bradley.  '•The  Fundamental  Solutions  fora 
Shallow  Shell  with  an  Arbitrary  Quadratic  Midsurface.  “ J.  Appl.  Mech..  43 
(Trans.  ASME.  92,  Series  E)  (1976).  p.  286. 


These  relationships  will  be  useful  in  obtaining  numerical  solutions  to 
the  problem. 

It  should  be  noted  that  the  bending  moments  and  M^,  which  are 

dependent  on  the  changes  in  curvatures  d/?/dr  and  (fif  r),  are  always  singular 

under  the  force  F,  but  the  membrane  forces  and  are  singular  only 

when  the  quantity  (P,  - P^)C  is  Konisero.  When  the  Poisson's  ratios  of  the 

two  metal  layers  are  equal  (i.e,,  Pj  - Plt  as  assumed  in  the  latter  part  of 

Ref.  1),  the  quantity  (J>,  - P^)  in  Eq.  (10)  above  vanishes,  and  the  stress 

resultant  N becomes  finite  under  the  force  F,  as  in  a single* layered  shallow 
r 7 

spherical  shell.  But  if  Pj  and  P^  are  slightly  different  (say,  P^  * 0.  HI  and 

P^  = 0.29,  as  used  in  the  numerical  examples  below),  N contains  a singularity 

•u  r 

under  the  force  F.  This  phenomenon  is  somewhat  similar  to  the  effect  of  a 
slight  deviation  of  the  shell  geometry  from  perfectly  spherical  (i.e,,  slightly 
ellipsoidal);  in  this  case,  a singular  term  proportional  to  the  difference  of  the 

,,  7 

two  principal  curvatures  exists  in  Nr  as  shown  by  Flugge  and  Filing. 

This  rather  strong  effect  of  the  condition  Pj  ■ P,  on  the  membrane  - 
stress  singularity  raises  a question  of  both  theoretical  and  practical  signifi- 
cance: Are  the  conclusions  drawn  from  Wittrick's  analysis  * (in  which  the 
above  condition  was  used)  concerning  the  thermoelastic  stability  of  the  bi- 
metallic curved  disk  strongly  dependent  on  that  condition'*  It  is  shown  below 
that,  when  the  condition  Pj  3 P,  is  relaxed,  all  the  major  conclusions  arrived 
at  in  Ref.  1 are  still  valid,  except  that  the  symmetry  property  of  the  solutions 
| see  Eq.  (**)  of  Ref.  l]  no  longer  holds  exactly. 

B.  FREE  DISK  INITIALLY  IN  THE  FORM  OF  A 
ggATlgW  SPHERICAL  CAP 

When  the  initial  form  of  the  disk  is  a shallow  spherical  shell,  its  initial 
reference  surface  may  be  approximated  by 


2 

This  gives  8 ■ 2h  r/a  . In  the  absence  of  the  central  force  F,  it  may  be 
o o 

readily  shown  by  direct  substitution  that  the  following  three  particular  solu- 
tions satisfy  Eqs.  (6)  and  (7)  and  all  boundary  conditions  exactly. 

Solution  (a).  Stress-free  State: 

0 ~ 0 and  N =0,  when  S - 0 
r 


Solution  (b),  Flat  State: 

0 - - 0o  and  Nr  « N rm(*  ’4). 
Solution  (c).  Reversed  State: 

P - ■ 2P0  and  Nr  * 0,  when 


when  B - B 

m 

e * is 

m 


In  the  above  solutions.  8 ^ is  the  "mean  temperature"  discussed  in  Ref.  1, 
and  may  be  obtained  by  substituting  0 « - 0q  in  Eq.  (9).  Thus 


« 


(1  * 

ZF 


(13) 


The  term  N is  the  maximum  compressive  stress  - resultant  at  the  disk 
center  when  the  disk  is  flat  (cf,.  Eq.  (34)  of  Ref.  l];  it  may  be  obtained  by 
substituting  0*  - 0 into  Eq.  (t>),  integrating  the  resulting  equation  and 
determining  the  integration  constant  by  Eq.  (8).  Thus 


N 


rm 


e 


(14) 


The  above  three  particular  solutions  for  a free  disk  are  essentially  the  same 
as  those  discussed  in  detail  in  Kef.  1,  except  that  Eqs.  (13)  and  (14)  are  now 
also  valid  for  the  more  general  case,  When  Vj  and  are  set  equal, 

Eq.  (14)  agrees  with  Eq.  (34)  of  Ref.  I,  but  an  attempt  to  reduce  Eq,  (13) 
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above  to  Eq.  (46)  of  Ref,  1 will  reveal  that  the  coefficient  K of  the  latter 
equation,  at  given  by  Eq.  (47)  thereof,  contains  an  extraneous  factor  (1  4 V)t 
which  is  believed  to  be  an  algebraic  error.  This  error  invalidates  the  cal- 
culated values  of  K in  Table  3 of  Ref.  1,  but  does  not  affect  the  analysis  of 
Refs.  1 through  3 in  any  other  way. 

Consider  next  the  symmetry  property  of  the  solutions  with  respect  to 
the  particular  solution  (b)  above  (i.e.,  the  fiat  position  of  the  disk)  as  dis- 
cussed in  Eq.  (33)  of  Ref.  1.  Specifically,  if 


S - • $ +5  and  N = N , when  S * 63  - 6? 

^ r r r m 


is  a solution  to  the  free  disk  problem,  examination  is  required  as  to  whether 
the  following  is  also  a solution: 

B * - B - B end  N • T5  . when  9 • + F . (15) 

r T T m 


Direct  substitution  of  Eq.  (IS)  into  Eq.  (6)  indicates  that  Eq.  (6)  is  not  satis- 
fied, owing  to  the  presence  of  its  second  term.  Therefore,  this  symmetry 
property  of  the  solutions  no  longer  holds  exactly  if  pj  * l^,  although  It  may  be 
expected  to  hold  approximately  if  P|  and  are  close. 


Finally,  it  is  necessary  to  reexamine  the  critical  initial  height  of  the 

disk  h , which  is  the  minimum  value  of  h for  the  thermal  snap  action  to 
c r o j 

occur.  A detailed  examination  of  Witt  rick's  procedure*  for  determining  h 


cr 


shows  that  the  procedure  will  remain  valid  for  the  more  general  case  of 
unequal  Poisson's  ratio  if  the  following  generalised  definitions  of  the  para- 
meters X and  4^  are  used  in  place  of  Eqs.  (28)  and  (30)  of  lief,  1: 

(» - 4)* 1 1/2 


•2q» 

n + *,d>d* 


7T5T 


2h 


(‘  * *'b)B 
2TT* 


1/2 


(16) 


(17) 
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Note  that  the  Poisson's  ratio  in  Eqs.  (31),  (38)  and  (42)  of  Ref.  1 must  also  be 
changed  to  as  seen  from  Eq.  (23)  thereof.  With  these  minor  changes, 
the  critical  values  of  may  be  calculated  for  a given  value  of  P^  through  the 
same  procedure  described  in  Ref.  1.  This  calculation  has  been  performed 
for  various  values  of  P^,  and  the  results  are  listed  in  Table  1 for  reference. 


Table  1.  Critical  Values  of 


"d 

0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

*cr 

8.  1590 

8.4186 

8.6537 

8.8674 

9.  0627 

9.2417 

If  an  "equivalent  thickness"  of  the  disk  is  defined  as 

t • (12  D*,'B)1/Z 

eq 


(18) 


then,  from  Eq.  (17),  it  is  seen  that  thermal  snap  action  of  the  disk  will 
occur  if 


h«> 


0 t 
*cr  eq 


- *) 


(19) 


Condition  (19)  above  reduces  to  Condition  (44)  of  Ref.  1 when  P|  ■ 

From  the  above  discussion,  it  may  be  concluded  that  when  the  curved 
disk  is  made  of  two  material  layers  with  unequal  Poisson's  ratios,  its  thermo- 
elastic snap  behavior  remains  essentially  the  same  as  described  in  Ref.  1. 

When  a central  force  is  present  on  the  disk,  however,  no  simple, 
closed-form  particular  solutions  can  be  found  similar  to  those  found  for  a 
free  disk.  Furthermore,  the  problem  is  no  longer  tractable  by  following  the 
power- series  solution  procedure  success  hilly  employed  in  Ref.  2.  Therefore, 
a direct  numerical  solution  was  pursued  in  the  following  analysis. 
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IV.  NUMERICAL  SOLUTIONS 


The  common  numerical  technique*  for  solving  two-point  boundary- value 

q 

problem*  can  be  classified  as  either  " finite -difference"  or  "shooting" 
methods.  For  a nonlinear  problem,  such  as  the  present  one,  finite -difference 
methods  will  result  in  a large  system  of  nonlinear  algebraic  equations,  thus 
losing  much  of  their  attractiveness  because  of  the  difficulty  in  itself  of  solving 
a large  nonlinear  algebraic  system  on  a digital  computer.  The  shooting 
method  is  a procedure  in  which  the  differential  equations  are  numerically 
integrated  as  an  initial  value  problem  with  estimated  unknown  initial  values, 
and  the  boundary  conditions  at  the  far  end  are  satisfied  by  iterative  correc- 
tions of  these  initial  estimates.  In  the  following  discussion,  numerical  solu- 
tions will  be  obtained  with  t'..e  aid  of  an  existing  FORTRAN  IV  code  developed  ¥ 

by  Gold10,  which  is  based  on  an  improved  shooting  technique  due  to  Morrison, 

Q 

Riley  and  Zancanaro. 

In  order  to  use  the  shooting  method,  the  two-point  boundary  value  prob- 
lem is  restated  in  the  following  form  as  a system  of  four  first-order  differ- 
ential equations 


du 

TF  ' 


i 

r 


(20a  1 


(20b) 


(20c) 


(20d) 


“3 

Morrison,  D.  D..  J.P.  Riley.  andJ.F.  Zancanaro,  "Multiple  Shooting 
Method  for  Two- Point  Boundary  Value  Problems,  " Communications  of  the 
Assoc,  for  Computing  Machinery,  5^(1 0e2).  p.  613. 

IOGold,  Barbara  J..  "Subroutine  TPRV.  " The  Aerospace  Corporation, 
Mathematics  and  Computation  Center  Library  (l*t>8). 
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where  tC^  it  the  change  in  meridional  curvature,  <r  it  the  meridional  mem- 
brane etrain,  and  N^,  rctpectively,  are  the  ttreet-  and  couple- retultantt 
in  the  circumferential  direction.  Thete  four  auxiliary  variable!  are  related 
to  the  four  primary  variable!  fi,  u.  and  by  the  following  algebraic 
equation!.* 

=^[Mr-|NrtQ*e]+,,Df 

<r  =®[Nr-CK  *'cf)+pe]-,'B7 

(21) 

Nfl  ’ B(?  + VB*r)  ‘ c(f  * VC*r)  ■ P® 

M8  * * N6  - D*  (f  * *'D*r)  - Q*« 

The  boundary  condition!  at  the  ditk  center  are 

/?((>)  » 0;  u(0)  ■ 0.  (22) 

Thote  at  the  edge,  r * a.  are 

Nr(a)  « 0;  Mf(a)  * 0.  (25) 

Note  that  a conatant  force  F ■ - KqI  which  appear!  in  Eq.  (20d).  hat  been 
attumed.  Method!  of  treating  a variable  force  will  be  deacribed  in  Section  V-A. 

In  order  to  circumvent  the  numerical  problem  aaeoctated  with  the  aingu- 

larity  at  r « 0,  it  it  aa  turned  that  the  atrtker-pin  ha  a a a mall  but  finite  radiua 

r 5 R , , and  that  the  force  F ia  produced  by  a uniformly  distributed  prea- 
pin  o 

aure  within  the  amall  region  r S Since  the  a utiona  are  now  finite,  the 

condition  of  axial  aymmetry  give# 

N « N . and  M » (for  r-e0).  (24! 

" 0 r V 
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By  combining  the  above  conditions  with  Eqs.  (21),  it  can  be  easily  shown  that 


^ * - *r,  and  ^ (for  r-*0).  (25) 

With  the  aid  of  Eqs.  (24)  and  (25),  the  singularity  in  Eqs.  (20)  and  (21)  can  be 
removed.  Thus,  F.qs.  (20c)  and  (20d)  are  replaced  by 


and  Eqs.  (21)  are  replaced  by 

«r g rrr^~[Mr  • yNf  * Q*e] 

N.  « (I  ♦ tOB  « ♦ (1  ♦ r_tC*  ♦ P 6 

9 or  ^ r 

M0  ■ SNe  *(l  ♦VD*VC*« 


(for  r-*0) 


(for  r-0) 


(20c1) 


(20d‘) 
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From  extensive  numerical  results,  it  was  found  that  as  long  as  R « a. 

pin 

the  solutions  are  independent  of  the  assumed  values  of  R^—  as  for  as  the 
disk  center  height  h and  the  snapping  temperatures  0u  and  B j ore  concerned. 
However,  the  local  bending  moments  and  Mq  were  found  to  increase  slowly 


with  d< 


"pin  v*'  '••• 


A 


ESTIMATES  OF  THE  UNKNOWN  INITIAL  VALUES 


The  key  factor  in  an  efficient  application  of  the  shooting  method  is  to 
provide  reasonably  close  estimates  of  the  unknown  initial  values,  in  this 
case  Nr(0)  and  M^O),  If  a given  temperature  0^  is  sufficiently  below  the 
snapping  temperatures  [say,  0,  £ (1/2)6^,  the  following  estimates  are  close 
enough  to  yield  solutions  (accuracy  requirement:  four  significant  digits! 
within  eight  iterations. 

Nr(0t  * Nrm  (l  - - <*{.  - •'b  K-  U6» 

Mr,0,  „ — yP-5  Nrm  (l  - /T^J  - <1  ♦ <n(Rpin/*'  (27> 

Note  that  the  first  terms  of  Eqs.  (26!  and  (27)  are  obtained  empirically  by 
assuming  that  the  st res* - resultant  Nf  increases  with  the  temperature  B in  a 
parabolic  form,  while  the  second  term*  thereof  were  obtained  from  the 
singular  solutions,  Eqs.  (10)  and  (11), 

When  final  solutions  at  two  selected  temperatures  0^  and  0 ^ have  been 
obtained  using  the  above  initial  estimates,  subsequent  solutions  for  higher 
temperatures  may  start  with  better  estimates  of  the  unknown  initial  values  by 
using  linear  extrapolation.  A sequence  of  solutions  may  then  be  obtained  at  a 
preselected  temperature  increment  40  until  no  convergence  can  be  reached. 
This  will  indicate  that  the  upper  snapping  temperature  0 > has  baen  exceeded. 

At  this  point,  40  is  reduced  by  a factor  of  ten.  and  the  solution  is  carried 
forward  from  the  last  convergent  temperature.  This  procedure  is  repeated 
until  the  upper  snapping  temperature  0u  is  determined  to  the  desired  accuracy. 
The  lower  snapping  temperature  0*  is  determined  in  a similar  manner  with 
a temperature  decrement  -40.  A computer  code  has  been  developed  to  carry 
out  the  above  procedures  automatically. 


B.  STRAIN  ENERGY  CALCULATIONS 


In  order  to  investigate  the  stability  of  the  various  equilibrium  positions, 
the  strain  energy  in  the  disk  is  calculated  lor  each  convergent  solution.  Let 
<U  denote  the  strain  energy  density  per  unit  area  ol  the  reference  surface; 

then 

«U  - 1 j <Vr  +Ve)df  <28) 

■h 

where  ( is  the  thickness  coordinate  measured  from  the  bond  surface  of  the 
two  layers.  The  elastic  strains  e^.  and  e^  are  given  by  the  relationships 

e'  ■ e -0.6  ■ («  + S*t 

(29) 

•g  ’ *0  -°i*  * (€0  -aia  + 


Note  that  assumes  the  value  Oj  for  ( > 0,  and  Ctj  *or  t < quantities 

< • u/r  and  #0  e -fi/r  are,  respectively,  the  membrane  strain  and  curva- 

ture change  in  the  hoop  direction.  Substituting  Eqs.  (29)  Into  (28)  and 
integrating  the  result,  one  obtains 
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Finally,  the  stresses  a and  (7*  can  be  expressed  in  terms  of  the  membrane 

r v 

strains  and  €q  and  curvature  changes  Kr  and  Kg[cf.,  Eqs.  (4)  of  Ref.  l]. 
Equation  (30)  can  thus  be  simplified  to  read  as  follows: 


au  -i[Nr«r  + Ve  +MA  + Ve1 

+ p2e2  - ±peur  +<0)  - |oe(*r  + *6>  (3»» 


where 


e,4i  e242 
2 = 1 -•'»  ‘-"z 


(32) 


Integration  of  the  above  expression  for  dU  over  the  entire  reference  surface 
gives  the  total  strain  energy  in  the  disk 


a 2 7T 

"<  •//  • 

0 ^0 


6 U rd0dr  = 2?r/  6 U rdr 

■'o 


4 


(33) 


Note  that  the  work  done  by  the  force  F,  which  gives  the  strain  energy  increase 
in  the  armature  from  the  preload  state,  must  be  added  to  to  obtain  the 
total  strain  energy  of  the  disk-armature  system.  That  is 

he 

U,o«a.  * Ud  7 ' F(h)  dh  <34) 

hl 

where  h.  and  h are  the  disk  center  height  at  the  striker- pin  engagement 
i e 

position  and  at  the  equilibrium  state  under  consideration,  respectively. 
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For  a given  equilibrium  state  at  a specified  temperature,  the  total 
strain  energy  calculated  from  Eq.  (34)  must  be  a stationary  value  with  respect 
to  a virtual  displacement  dh  (corresponding  to  an  admissible  change  of  disk 
shape  dw);  it  may  be  a minimum,  a maximum,  or  a point  of  inflection.  The 
equilibrium  state  is  stable  only  when  the  total  strain  energy  is  a local  mini- 
mum. If  the  strain  energy  is  a local  maximum,  the  disk  is  in  unstable 
equilibrium,  and  it  will  snap  to  one  of  two  possible  stable  equilibrium  posi- 
tions, having  lower  total  strain  energy.  If  the  strain  energy  is  at  an  inflec- 
tion point,  the  disk  is  again  at  unstable  equilibrium,  and  it  will  snap  to  a 
unique  stable  equilibrium  position  (see  Fig.  5 for  additional  clarification  of 
these  three  types  of  equilibrium  positions). 


t 
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V.  NORMAL  AND  ANOMALOUS  SWITCHING  CYCLES 


In  order  to  investigate  the  switching  behavior  of  the  type  (a)  BMD 
thermostats,  a detailed  numerical  analysis  was  performed  on  a typical  BMD 
thermostat  design  which  has  been  used  extensively  in  recent  space  programs 
and  has  experienced  dithering  problems.  Its  parameters,  physical  and  cal- 
culated, and  certain  performance  specifications  are  shown  in  Table  2. 

Although  the  nominal  switching  band  of  the  particular  batch  of  thermostats 
being  analyzed  is  listed  as  5.6°C  (10°F),  it  is  interesting  to  note  that  the 
specification  gives  the  switching  band  as  1.  l°C-9.  0°C  (2°F-16°F),  an  indi- 
cation of  wide  unit-to-unit  variations. 

As  discussed  in  Ref.  3,  the  bimetallic  disks  used  in  these  thermostats 
are  formed  by  a series  of  stamping  operations,  with  only  incomplete  stress 
relief;  consequently,  the  end  product  undoubtedly  contains  high  residual 
stresses  at  room  temperature.  One  evidence  of  the  presence  of  these  residual 
stresses  in  the  disks  is  the  subcritical  disk  center  height  at  room  temperature. 
The  stress-free  state*,  even  if  one  assumes  it  can  be  recovered  by  a temperature 
change  alone,  is  at  an  unknown  temperature,  which  cannot  be  measured  directly. 

In  the  following  applications  of  the  theory,  the  "inverse"  procedure  first 
used  in  Ref.  3 is  employed  to  determine  the  initial  disk  center  height  hQ;  i.  e., 
a series  of  trial  values  of  hQ  is  used  in  the  analysis  until  the  calculated 
snapping  temperature  band  Ad  ■ 0u  - 6 1 matches  the  measured  band  for  a 
free  disk,  AT  ■ Ty  - Tj  (note  that  the  latter  is,  in  general,  different  from 
the  thermostat  switching  band  AT',  as  will  be  discussed  in  more  detail  in  the 
next  section). 

The  rationale  for  this  procedure  is  based  on  the  fact  that  the  residual 
stresses  in  the  bimetallic  disk  must  satisfy  the  same  set  of  shell  equations 
and  boundary  conditions  as  do  the  thermal  stresses  caused  by  a temperature 
rise.  Whereas  the  thermal  stresses  are  caused  by  unequal  thermal  expan- 
sions in  the  two  layers  fi,e„  the  PS  and  Q+0  terms  in  Eqs,  (9)  through  (12) 
of  Ref,  1],  the  residual  stresses  are  primarily  caused  by  the  unloading  of  the 


Table  2.  Thermostat  Parameters 


a.  Disk  Physical  Parameters 

Thermal  expansion  coefficient 

“i 

* 26.3  X 10"6  0C_1  (14.6  X 10‘6  “F"1) 

“2 

» 2.2  X Kf6  “C"1  (1.2  X 10"6  °F’1) 

Young's  Modulus 

Ei 

■ 124.  1 GPa  (18.0  X 106  psi) 

E2 

a 156.5  GPa  (22.7  X 106  psi) 

Poisson's  Ratio 

"l 

= 0.30 

"2 

= 0.29 

Thickness 

a 0.  12  mm  (0.0047  in.  ) 

*2. 

< 0.  11  mm  (0.0043  in. ) 

Disk  Radius 

a ■ 

6.35  mm  (0.25  in.) 

Disk  center  height  at  RT 

WRT  ■ mm  (0*  012  in.) 

b.  Performance  Specifications 

Upper  switching  temperature 

T 

u 

• 51.7  4 1.  7°C  (125  4 3°F) 

Lower  switching  temperature 

ri 

• 46.1  °c  <*15  °F) 

Switching  band 

AT 

■ 1.  l°C  - 9.  0°C  (2°F  - 16°F) 

Striker-pin  force  required  to 
overcome  contact  preload 

FP 

- 0.68  N 

c.  Calculated  Parameters 

Critical  value  of  h 

o 

h 

cr 

» 0.435  mm  (t*. 01714  in.) 

Disk  center  height  at  "stress- 
free"  state 

N. 

• 0.48  mm  (0.0189  in.) 

Mean  temperature  of  disk  snap 
points 

• 150.8  (°C  from  reference  temperaturel 

Upper  snap  point 

®u 

• 156.2  (°C  from  reference  temperature) 

Equivalent  Poisson's  ratios 

"b 

* 0.295.  p(.  • 0.0912.  PD  • 0.295 

Spring  constants  of  armature  at 
striker -pin* 

kl 

■ 26  N/mm.  k^  * 5 N/mm 

Calculated  from  a finite -element  plat*  analysis  and  verified  by 
lab  measurements 
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unequal  yield  stresses  in  the  two  metal  layers  at  the  end  of  the  forming 
process.  Therefore,  within  the  scope  of  a thin  shallow  shell  theory,  it  is 
not  unreasonable  to  expect  that  this  residual  stress  state  may  be  approxi- 
mated by  a thermal  stress  state  to  the  extent  that  their  effects  on  such  gross 
structural  behavior  as  snap  buckling  or  free  vibration  would  be  equivalent. 

Application  of  the  above  procedure  indicates  that,  for  the  thermostat 
discussed  herein,  the  disk  snapping  band  of  SB  - 10,B°C'  corresponds  to  the 
initial  center  height  hQ  * 0.48  mm,  and  that  the  "equivalent  stress-free  state" 
is  at  the  reference  temperature  T * - 104,  5°(*. 

A . NORMAL  THERMOSTAT 

The  effects  of  a central  force  on  the  thermoelastic  deformation  of  the 
disk  is  shown  in  Fig.  4,  in  which  a family  of  curves,  corresponding  to  a set 

of  constant  forces  (Fq  = 0,8  F^,  F^ 2.5  F^)  is  plotted,  together  with 

the  baseline  curve  for  the  free  disk  (Fq  = 0),  These  S- shaped  curves  show 
the  disk  center  height  h of  the  equilibrium  configuration  as  a function  of  the 
disk  temperature  B . An  interesting  observation  from  Fig.  4 is  the  striking 
similarity  of  all  these  curves.  Within  the  accuracy  of  the  numerical  results, 
it  was  found  that  each  of  these  curves  is  congruent  to  the  baseline  curve  and 
is  at  a temperature  shift  SB  from  this  curve.  The  shift  SB  was  found 
empirically  to  be  proportional  to  the  magnitude  of  the  force  F ; i.e.. 


$0  « 0.  HO  (F  /Q*L  OM 

o 

The  above  result  indicates  that  the  snap  band  SB,  the  two  values  of  disk 
center  height  at  which  instability  occurs,  and  the  throw  of  the  disk  during 
snap-through  are  not  affected  by  the  presence  of  the  stresses  due  to  a constant 
central  force  Fq,  This  feature  provides  farther  support  to  the  inverse  pro- 
cedure used  previously  in  the  treatment  of  the  unknown  residual  stress  con- 
dition, since  it  shows  that  the  snapping  characteristics  of  the  disk  are  not 
sensitive  to  the  detailed  preetress  distribution  therein. 


Fig.  4.  Disk  Center  Height  versus  Temperature  — Normal 
Switch  (switching  cycle  ABC'D'G1) 
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The  above  family  of  h-0  curves  (for  different  FQ  values),  combined 
with  a family  of  vertical  straight  lines,  provides  a network  of  curvilinear 
coordinates  for  tracing  equilibrium  states  of  a disk  with  a variable  striker- 
pin  force  F = F(h).  If  the  striker-pin  is  assumed  to  engage  the  disk  at  its 
flat  position  (i.e.,  h^  = 0),  and  the  force  function  is  bilinear  as  illustrated  in 
Fig.  2,  the  h-0  curve  corresponds  to  the  heavy  line  ABUD'C'  in  Fig.  4. 

Note  that  this  curve  no  longer  has  a simple  S- shape,  but  has  abrupt  slope 
changes  at  both  the  point  of  striker-pin  engagement  1 and  that  of  contact 
reparation  J.  It  is  important  to  point  out  that  the  loop  G'BC'D'G’  represents 
a normal  thermostat  switching  cycle.  In  this  cycle,  the  upper  snap  point  B 
remains  the  same  as  that  of  the  free  disk,  but  the  lower  snap  point  D1  is  dif- 
ferent from  point  D of  the  free  disk.  It  is  apparent  that  the  lower  snap 
temperature  may  be  raised  significantly  by  the  interaction  force  F;  conse- 
quently, the  snap  band  for  the  switch  AS*  may  be  considerably  less  than  AS , 
the  free -disk  snap  band.  However,  the  disk  center  height  at  which  the  lower 
snap  action  occurs  appears  to  be  only  slightly  affected  by  the  force  F. 

The  total  strain  energy  of  the  disk-armature  system,  as  well  as  that  of 
the  free  disk,  are  illustrated  in  Fig.  5 for  selected  temperatures.  I;  is 
worthy  of  note  that  in  the  neighborhood  of  the  lower  snap  point  D',  the  total 
strain  energy  decreases  only  slightly  as  the  disk  center  moves  toward  the 
contact  closing  point  Gp  an  indication  that  the  lower  snap  action  of  this 
thermostat  may  start  with  an  undesirable,  slow  motion:  such  motion  could 
give  rise  to  hesitant  contact  closure. 

B.  ANOMALOUS  THERMOSTATS 

It  is  apparent  from  an  inspection  of  Fig.  4 that  the  position  of  striker - 
pin  engagement,  or  h^,  must  be  kept  within  a very  tight  tolerance  for  the 
thermostat  to  function  properly.  For  the  thermostat  under  consideration, 
this  allowable  tolerance  is  on  the  order  of  0. 1 mm  (0.004  in. ).  Since  control 
of  manufacturing  tolerances  to  this  level  is  difficult,  if  not  impractical,  it  is 
not  surprising  that  high  rejection  rates  in  acceptance  tests  and  persistent 
problems  during  service  were  experienced  with  this  thermostat. 
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Fig.  5.  Strain  Energy  ir  Free  Ditk  (dashed  Unes)  and  in 
Disk- Armature  System  (solid  lines)  — Normal 
Switch  (see  Fig.  4) 
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First,  consider  the  thermostat  in  which  the  striker-pin  is  too  long, 
so  that  at  the  upper  snap  temperature  it  engages  the  disk  before  the  latter 
snaps  past  the  flat  position.  In  such  a case,  the  drop  in  strain  energy  from 
the  upper  snap  point  B to  the  contact -opening  point  Cg  (see  Fig.  5)  is  rela- 
tively small.  This  available  energy  may  be  insufficient  to  overcome  the 
barrier  of  electric  potential,  and  the  contacts  may  fail  to  open  immediately. 
Although  further  temperature  rise  will  increase  the  disk-pin  interaction  force 
and  eventually  open  the  contacts,  this  delay  in  contact  opening  is  undesirable 
because  it  will  leav<  *>e  contacts  momentarily  in  a condition  of  poor  contact, 
thus  susceptible  to  arcing. 

Consider  next  the  thermostat  in  which  the  striker-pin  is  too  short,  so 
that  at  the  lower  snap  temperature  the  contacts  close  before  the  disk  is  set  in 
motion  (*ee  Fig.  6).  In  such  a case,  the  disk-armature  system  is  in  a stable 
condition  when  the  contacts  are  in  a threshold  contact;  consequently,  the 
contacts  will  alternatively  open  and  close  with  thermally  induced  motion 
instead  of  snap  action.  This  anomalous  form  of  switching  will  be  referred  to 
in  this  report  as  "creepage". 

It  should  be  pointed  out  that  the  above  anomaly  may  escape  detection  in 
a functional  test  in  which  the  temperature  changes  are  controlled  independently 
of  the  electrical  contacts  (such  as  by  a liquid  bath,  a commonly  used  method), 
since  the  disk  temperature  0 would  then  continue  to  decrease  after  contact 
closing.  When  the  temperature  reaches  point  1 (Fig.  6),  the  disk  disengages 
from  the  striker-pin  and  snaps  to  the  prebuckling  configuration  G\  However, 
if  a sufficiently  slow  rate  of  change  of  the  disk  temperature  is  maintained  in 
such  a test,  a time  lapse  should  be  easily  detected  between  the  contact  closing 
(which  may  be  recorded  as  an  electric  signal)  and  the  disk  snapping  (an 
acoustic  or  accelerometer  signal). 

If  the  thermostat  is  used  to  control  power  to  a heater  upon  which  the 
thermostat  is  mounted,  and  its  disk-pin  engagement  point  I is  made  slightly 
closer  to  the  base-plane  of  the  disk  than  in  the  previous  case  (as  shown  in 
Fig.  7),  a more  complex  and  interesting  switching  anomaly  may  result. 
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The  lower  instability  point  D (see  Fig.  7)  is  now  within  the  temperature  range 
of  the  transition  segment  IJ  of  the  h-0  curve;  therefore,  when  the  disk  is 
cooled  from  point  C'  to  point  D1,  it  will  move  toward  a stable  "transition 
equilibrium  state"  H,  at  which  the  contacts  are  closed  but  the  contact  force  is 
less  than  the  set  preload.  If  sufficient  heat  is  transmitted  from  the  heater  to 
reverse  the  temperature  decline  before  point  I is  reached,  the  disk  tempera- 
ture will  rise  toward  point  J along  IHJ  (see  Fig.  7);  at  J the  disk  becomes 
unstable  and  moves  to  the  stable  equilibrium  state  K.  This  motion  switches 
the  heater  off  and  causes  the  disk  temperature  to  drop  toward  point  D\  Thus, 
the  thermostat  may  repeatedly  perform  the  secondary  switching  cycle 
KD'HJK,  which  keeps  the  disk  temperature  "dithering"  within  a narrow  range 
(about  0.6°C,  as  shown  in  Fig.  7)  without  letting  the  disk  snap  to  a pre  buck  ling 
configuration. 

The  above  secondary  switching  loop  cannot  be  completed  if  the  disk 
temperature  is  controlled  externally  or  if  there  is  a sufficient  time  lag  be- 
tween contact  closing  and  disk  heating.  In  such  cases,  the  disk  temperature 
continues  to  drop  along  JHI  until  it  reaches  point  I.  At  this  point,  the  disk 
snaps  to  the  pre  buck  ling  configuration  C'and  completes  a primary  (or  designed) 
switching  cycle  G'BC'D'HIG*. 

Figure  7 illustrates  that,  at  a temperature  within  the  dithering  range 
(i.e.,  149. 14  < 0 < 149.  73),  the  disk-armature  system  has  five  equilibrium 
states.  In  order  to  investigate  the  stability  of  these  equilibrium  states,  the 
values  of  the  total  strain  energy  at  these  states  were  calculated  and  are  shown 
in  Fig.  8.  It  was  found  that,  during  the  secondary  switching  motions,  the  total 
strain  energy  varies  by  such  a slight  amount  that  an  additional  significant 
figure  of  accuracy  was  needed  in  the  numerical  calculations  to  obtain  reason- 
able resolution.  This  indicates  that  during  the  secondary  switching  cycle, 
the  disk-armature  system  is  in  a state  of  nearly  neutral  equilibrium  if  only 
the  thermoelastic  forces  are  present.  Consequently,  any  effect  that  would  be 
normally  inconsequential  in  a stable -equilibrium  region  could  assume  an  im- 
portant role  in  the  secondary  switching  cycle.  Such  factors  as  inertial  forces 
due  to  shock  or  vibration,  a slight  asymmetry  of  the  oisk,  material  aniso- 
tropy, or  electrostatic  force  due  to  voltage  across  the  sontacts  could  signifi- 
cantly alter  the  shape  of  the  secondary  switching  cycle  KD'HJK  (see  Fig.  7). 
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Fig.  0.  Strain  £atrgy  in  Frn  Disk  (dashtd  lints)  and  in 
Disk-Armaturt  Systtm  (solid  lints)— Anomalous 
Switch  (stt  Fig.  ?) 
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VI.  STRESSES  IN  BIMETALLIC  DISK 


It  was  shown  in  Ref.  S that  high  mean  and  cycli . stresses  may  develop 
in  a typical  bimetallic  disk  during  its  thermal  snar  a*  , eras.  Therefore,  it  is 
desirable  to  design  a BMD  thermostat  with  a view  to  tmiting  the  maximum 
mean  and  cyclic  stresses.  Since  Ref.  3 did  not  include  the  effects  of  the 
disk-pin  interaction,  a detailed  investigation  was  performed  in  the  present 
analysis  to  study  whether  this  interaction  force  increases  or  decreases  the 
severity  of  these  stresses. 


For  normal  thermostat  behavior,  illustrated  in  Figs.  4 and  5,  the 

maximum  cyclic  str for  both  the  free  disk  and  the  disk  with  pin-force 

are  shown  in  Fig.  9 a.  functions  of  the  disk  temperature  9.  These  maximum 
cyclic  stresses  occur  on  the  surfaces  of  the  disk,  at  its  center.  It  is  signifi- 
cant to  observe  that  the  striker-pin  force  reduces  the  cyclic  stress  ranges  of 
both  material  layers,  evidently  because  it  reduces  the  throw  of  the  snap 
actions.  Therefore,  from  the  viewpoint  of  material  fatigue  damage,  the 
presence  of  the  disk-pin  interaction  force  does  not  appear  to  increase  the^ 
severity  of  the  cyclic-stress  problem  from  that  of  the  baseline  free  disk. 

It  should  be  pointed  out.  however,  that  the  above  result,  are  based  on  the 
assumption  that  the  disk-pin  interaction  force  is  distributed  over  a finite  area 
with  a radius  approximately  equal  to  the  disk  thickness;  in  addition,  the 
impact  or  dynamic  effects  of  the  snap  action  were  neglected. 

The  maximum  tensile  stresses,  which  are  usually  the  more  critical 
factor  in  determining  crack  development  in  the  material,  are  shown  in  Fig.  10 
for  the  same  disk,  with  and  without  the  striker-pin  force.  These  maximum 
tensile  stresses  occur  at  the  edge  of  the  disk  a.  hoop  stresses  and  always  in 
material  layer  2 (the  low  expansivity  metal).  Figure  10  indicate,  that  the 
value  of  the  maximum  tensile  stress  is  not  affected  by  the  presence  of  the 
disk-pin  interaction,  since  it  occurs  just  as  the  disk  snaps  past  the  flat 
position  during  both  the  upper  and  lower  snap  actions. 


STRESS  (MPa)  STRESS  (MPa) 


(b)  Material  2 Surface  Stress  at  Disk  Center 

Fig.  9.  Comparison  of  Maximum  Cyclic  Stresses  in  Free  Disk 
(dashed  lines)  and  in  Disk  with  Pin- Force  (solid  lines) 
--Normal  Switch  (s?e  Fig.  4) 
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TEMPERATURE,  8 <°C> 


Fig.  10.  Comparison  of  Maximum  Tonsils  Strossos  in  Froo  Disk 
(dashod  linos)  and  in  Disk  with  Pin- Fores  (solid  linos)- - 
Normal  Switch  (boo  Fig.  4) 
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To  summarize  the  effects  of  the  striker-pin  force  on  the  disk  stresses, 
it  may  be  stated  that,  other  than  the  possible  transient  effects  due  to  disk-pin 
impact,  the  maximum  tensile  and  cyclic  stresses  in  the  disk  during  the 
switching  operation  of  a thermostat  are  bounded  by  those  in  the  free  disk 
during  its  snap  cycle. 
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VII.  CONSIDERATIONS  FOR  THERMOSTAT  DESIGN  IMPROVEMENT 


From  the  above  analysis,  it  may  be  seen  that  dithering  or  creepage  of 
BMD  thermostats  results  from  a combination  of  the  following  three  conditions: 

1.  A relatively  high  disk-pin  interaction  force  Fp  for  the  contact 
opening,  creating  a significant  region  of  '‘transition11  equilibrium 
(segment  IJ  of  the  curves  shown  in  Figs.  4,  6 and  7) 

2.  A relatively  short  striker-pin,  causing  its  disengagement  from 
the  disk  at  a position  (h  = hj)  that  is  close  to  the  lower  switching 
point 

3.  A slow  rate  of  temperature  change  of  the  environment  and  a rapid 
thermal  response  of  the  switch-heater  system,  so  the  disk 
temperature  does  not  experience  appreciable  overshoot  from  the 
switching  temperatures 

Among  the  above  conditions,  the  last  is  partly  imposed  by  the  thermal  environ- 
ment and  partly  by  a characteristic  of  the  thermal  system  in  which  the  ther- 
mostat is  used;  therefore,  it  cannot  be  eliminated  by  change  of  the  thermo- 
stat design  alone.  Although  it  has  recently  been  suggested  that  a delayed 
response,  and  thus  an  overshoot,  of  the  thermostat  temperature  be  purposely 
designed  into  the  thermal  system  to  eliminate  temperature  dithering,  it  is 
not  a desirable  solution  for  the  following  two  reasons.  First,  for  a given 
thermal  capacitance  between  the  heater  and  the  thermostat,  the  temperature 
overshoot  will  be  sufficient  to  prevent  dithering  only  above  a certain  rate  of 
temperature  change;  if  the  rate  is  slower,  dithering  will  occur.  Second, 
and  more  importantly,  this  technique  only  attempts  to  eliminate  the  "symptom" 
(the  temperature  dithering),  but  will  not  alleviate  the  underlying  problem, 
i.c.,  the  nearly  neutral  equilibrium  and  poor  contact  condition  at  the  lower 
switching  point.  Therefore,  improvement  of  the  thermostat  performance  and 
reliability  must  be  sought  from  design  changes  that  eliminate  the  first  two 
conditions  above. 
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Condition  1 may  b©  changed  either  by  reducing  the  contact  preload  or  by 
increasing  the  snap-temperature  band  AS  of  the  free  disk  (the  latter  does  not 
necessarily  imply  an  increased  switching  band  of  the  thermostat).  Since  a 
sufficient  contact  preload  is  necessary  tc*  provide  a good  electrical  contact 
and  to  prevent  chattering,  it  does  not  appear  promising  that  the  contact  pre- 
load in  existing  designs  can  be  significantly  reduced.  Consequently,  improve- 
ment can  best  be  achieved  by  using  bimetallic  disks  with  wider  snap- 
temperature  band.  From  Fig.  4 and  Eq.  05),  it  can  be  deduced  that  if 
AS  <0.  26  (F  /Q*)t  the  disk  will  always  yield  a bad  switch.  The  following 
condition  appears  to  provide  a suitable  criterion  for  the  selection  of  bimetallic 
disks: 


AS  £ 0.8  (F  /Q*).  (W 

r 

With  the  above  selected  values  of  a normal  switch  should  result  if  the  disk 
pin  engagement  point  l shown  in  Fig.  4 is  kept  within  the  middle  third  of  the 
temperature  range  AS* 

C ondition  2 above  may  be  alleviated  by  maintaining  a high  dimensional 
precision  for  the  parameter  hj  during  the  assembly  of  the  thermostats  or  by 
selecting  disks  with  increased  throw  of  the  snap  actions  so  that  the  admissible 
range  of  hj  is  significantly  increased.  In  order  to  ensure  that  the  actual 
value  of  h|  is  well  within  the  admissible  range,  detailed  analysis  should  be 
performed  for  each  thermostat  design,  and  tolerance  requirements  for  per- 
tinent dimensions  should  be  estimated  and  compared  to  practical  limits.  The 
following  formulas,  derived  from  a parametric  study,  enable  one  to  calculate 
the  initial  disk  center  height  h^  directly  from  the  measured  snap-temperature 
band  AS  without  going  through  the  t rial-and-e rror  process  mentioned  above. 

1-et  ( denote  a dimensionless  parameter  describing  the  initial  disk 
center  height,  i.e., 

< • r*  - ' " $r  - 1 <'7> 

cr  rcr 
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An  extensive  parametric  study  and  curve  fitting  yield  the  following  relation 
for  the  dimensionless  quantity  IT  - l) 

(T  - 1>  = l.07p/Z  (1  - 0.  2?t>{  - 0.087{Z)  (181 

where  T ~ & l&  , as  defined  in  Hef.  1.  Equation  (18)  is  valid  for  both  a 
u u ^m 

free  disk  and  a disk  with  a small,  constant  central  force,  provided  that  h 

o 

is  within  the  scope  of  a shallow  shell  theory  (say,  h < a/4)  and  that  the  two 
Poisson's  ratios  and  V , are  not  too  widely  different  (say,  | |/^  - |/,  | < 0,1), 
in  such  cases,  the  two  quantities  (l  - Tj  ) and  (r  - I)  may  be  considered  equal, 
and  the  snap-tempe rature  band  is  given  by  the  following  relationship: 
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Prom  Pq.  (ID,  it  may  be  readily  shown  that 
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Substitution  of  Kqs,  (18)  and  (40)  into  Pq.  (1‘M  yields 

(l  * V. .)l>*h  w » ■»  * 

AS  - -~4..'8;  it  * 0.7t»4{  0.  DM*-  O.OS74).  t4D 

Q a* 

When  the  measured  value  of  is  known,  Pq.  (41)  may  he  solved  for  the 

unknown  parameter  4.  which  may  then  be  used  in  Pqs.  (17)  and  (40)  to  cal 

culate  h and  S , respectively, 
o m r 

Finally,  it  is  important  to  remark  that  a better  understanding  and 
control  of  the  instability  of  the  relevant  material  properties  (e.g..  those  of 
INVAK  are  well  known  to  be  unstable)  and  alleviation  o(  the  residua)  stresses 
in  the  bimetallic  disk  must  accompany  any  design  improvement  effort;  other 
wise,  any  improvement  in  the  switch  performance  and  reliability  achieved 
through  design  changes  and  dimensional  precision  can  easily  be  negated  by  a 
drift  in  these  material  characteristics. 


VIII.  CONCLUSIONS 


The  analysis  and  numerical  examples  presented  in  this  report  provide  a 
more  detailed  understanding  of  the  operation  of  the  BMD  thermostat  by  includ- 
ing the  effects  of  the  disk-armature  interaction  force.  It  is  shown  that  the 
temperature  creepage  or  dithering  caused  by  the  anomalous  switching  actions 
of  a miniaturized  BMD  thermostat  may  be  attributed  to  the  improper  design 
of  this  interaction,  and  that  these  anomalies  may  escape  detection  in  an 
acceptance  test  if  the  rate  of  change  of  the  disk  temperature  is  not  sufficiently 
slow.  A few  theoretical  guidelines  have  been  suggested  for  design  improve- 
ment. 


SYMBOLS 


a 

b,c,d 

D* 


ErE2 

Vee 

F, 


P 

h ■ u( o) 
ho  * -%{0) 


Mr-Me 

Nr.Ne 

Q_ 


r 

T 

t 


*!•*! 


U« 

U 

U 

w 


total 


radius  of  the  edge  of  disk 

disk  constants  defined  in  Ref.  1 

D - C / B,  equivalent  bending  rigidity  of  disk 

Young’s  moduli  of  disk  materials 

total  strains  in  radial  and  circumferential  directions 
disk-striker-pin  interaction  force 

the  interaction  force  required  to  overcome  contact  preload 

disk  center  height  measured  from  base-plane 

initial  disk  center  height  (at  stress -free  state) 

disk  center  height  at  striker-pin  engagement  position 

•tress-couple  resultants  in  disk 

stress  resultants  in  disk 

transverse  shear  stress  resultant  in  disk 

radial  distance  from  disk  center 

temperature 

total  thickness  of  disk 

thickness  of  layers  1 and  2 

•train  energy  in  disk 

total  strain  energy  in  disk-armature  system 

radial  displacement  of  disk  reference  surface 

axial  distance  from  the  reference  surface  of  the 
disk  to  its  base -plane 
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SYMBOLS  (Continued) 


w 


w 


ar«2 

0 s dw  /dr 
© 

tr‘e$ 


t 

X 


axial  displacement  of  disk  reference  surface 
initial  shape  of  the  disk  (at  stress-free  state) 
coefficients  of  thermal  expansion  of  disk  materials 
rotation  of  disk  normal 

elastic  strains  in  radial  and  circumferential  directions  of 
disk  reference  surface 

thickness  coordinate,  measured  from  reference  (bond) 
surface  of  disk 

disk  constant  defined  in  Eq.  (16) 


V Poisson's  ratios  of  disk  materials 

t - (h  /h  ) - 1,  dimensionless  parameter  describing 

init&l  Curvature 


‘V'o 
t = e/e 


m 


stresses 

dimensionless  temperature  parameter 

variables  describing  the  deformed  and  initial  shapes  of 
disk 


e 


disk  temperature  measured  from  a reference  temperature 
at  its  stress-free  state 


upper  and  lower  snap  temperatures  of  disk 
0m  mean  temperature  of  disk  snap  points 


V 


changes  of  shell  curvature 
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APPENDIX 


SUMMARY  OF  WITTRICK'S  DERIVATION 


Consider  a bimetallic  disk  of  radius  a and  in  the  form  of  a shallow 

shell  as  shown  in  Fig.  ' <P-  HI.  The  ,ayers'  who‘e  thickne!‘8e',  ar* 
t|  and  t,.  whose  Young's  moduli  are  E,  and  E,.  and  whose  Poison's  ratios 

are  It,  a~„d  and  having  coefficients  of  linear  thermal  expansion  «,  and  a,, 
are  assumed" to  be  in  perfect  bond.  Wittrick's  basic  assumption  is  that  there 
exists  a reference  temperature.  T„.  at  which  the  disk  is  free  from  internal 
stresses  and  which  has  a slightly  curved  shape  represented  by  the  funct.on 
w ,r!  Other  assumptions  are  that  t uc  a and  that  the  magnitude  of  wp  is  the 
same  order  as  t.  (t  • t,  ♦ M-  Then,  based  on  the  large  deflection  theory  of 
thin  plates  and  shallow  shells,  the  radial  and  circumferential  strains  may  be 
expressed,  respectively,  as 


d 

d 


u 

r 


%.  dw 
S £ 

r dr 


(A- 1> 


where  u and  w are  the  thermoela.tic  displacement  component,  in  r-  and  »■ 
direction,  respectively,  of  a point  at  a distance  ( from  the  reference  (or  bond. 

surface,  and 


♦ w 


(A..M 


is  the  final  (or  deformed)  shape  of  the  disk. 

Removal  of  the  thermal  strain,  due  to  a temperature  change,  t?  - T - 1\,. 
give*  the  elastic  strains  due  to  stresses 


- e. 


Cl^. 


*0 


axe 


tA.M 


A-  l 


Note  that  a.  assumes  the  value  for  ( >0,  and  Ct^  for  ( < 0,  The  radial 
and  circumferential  stresses  are  then  given  by 


<rr  « V (e'+M'l 

r 1 - l/f  r 1 


°0  = 


i - 


(A  - 4 ) 


(A -5) 


V*e t N and  be  the  radial  and  circumferential  stress- resultants,  respec- 
r p 

tively 


V* 


(A-6) 


and  M and  M~  be  the  radial  and  circumferential  stress -couple  resultants, 
r tr 

respectively 

Mr  * I Or  dsk.  - J eei4Z  (A-; 

/».  /t. 


For  given  di»k,  the  following  quantities  were  defined: 


Ei  ■ 


l - V \ 


B • + E^t,, 


Ez  * 


B*  - *'iEiti + Wz- 


‘ (t>t\  - E‘tj|  C-  ■ ‘ ^EJtf  - F2E’tf). 


D s 1^1*1  + E2*z)’ 
Wi 

p = Try7  TTv{- 


°'  • i ('»*'*»  * V'4 

1 a^t2  \ 

Qsn'TrTi'^r7Tr)' 


(A-8) 


Then,  ueing  Eqe.  (A-l)  through  (A-5),  Eq.  (A 


■b)  may  be  integrated  to  give 
P0  - C 


d2we  1 

T T ■ c 7 dr 
dr 


_ ldwe 

p«-  C?-5T-C 


d2W 


(A -9) 


and  Eq.  (A -7*  may  be  integrated  to  give 


Mr  * I 


\(°t )»■(”-£)  ^r-(D' 


dw« 

"fls'sr 


CC*\  1 e 


/ 


(A-U) 


(A  - 12) 


Equations  (A-9)  and  (A- 10)  can  be  solved  for  u and  (du/dr).  Upon  elimina- 
tion  of  u.  one  obtaine 


(A-1JV 


In  the  absence  of  external  forces,  the  equations  of  equilibrium  are 


d_ 

dr 


57  <rPV  * N* 


(rMp)  - ♦ rN 


S 0 

dw 
r dr 


0 


(A  - 14) 

(A  - 15) 


Substitution  of  Eq.  (A- 14)  into  Eq.  (A- 13)  gives 


( A - 1 6 ) 


Also,  substituting  Eqs,  (A- 11)  and  (A- 12)  into  Eq.  (A- 15)  and  ueing  Eq.  (A- 14) 
gives 


(A- 17) 


Equations  (A  - lo I and  (A- 17)  are  two  coupled  differential  equations  governing 
Nr  am!  dw /dr. 


A-4 


